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Abstract 

The sensitivity achievable by a pair of VIRGO detectors to stochastic and isotropic grav- 
itational wave backgrounds produced in pre-big-bang models is discussed in view of the 
development of a second VIRGO interferometer. We describe a semi-analytical technique 
allowing to compute the signal-to-noise ratio for (monotonic or non-monotonic) logarithmic 
energy spectra of relic gravitons of arbitrary slope. We apply our results to the case of two 
correlated and coaligned VIRGO detectors and we compute their achievable sensitivities. 
We perform our calculations both for the usual case of minimal string cosmological scenario 
and in the case of a non-minimal scenario (originally suggested by Gasperini) where a long 
dilaton dominated phase is present prior to the onset of the ordinary radiation dominated 
phase. In this framework, we investigate possible improvements of the achievable sensitivi- 
ties by selective reduction of the thermal contributions (pendulum and pendulum's internal 
modes) to the noise power spectra of the detectors. Since a reduction of the shot noise does 
not increase significantly the expected sensitivity of a VIRGO pair (in spite of the relative 
spatial location of the two detectors) our findings support the experimental efforts directed 
towards a substantial reduction of thermal noise. 



1 The problem and its motivations 



Every variation of the background geometry produces graviton pairs which are stochastically 
distributed and whose logarithmic energy spectra represent a faithful snapshot of the (time) 
evolution of the curvature scale at very early times |L|. Indeed, one of the peculiar features 
of stochastic graviton backgrounds is that their energy spectra extend over a huge interval of 
(present) frequencies. Since gravitational interactions are much weaker than electromagnetic 
interactions they also decouple much earlier. Therefore the logarithmic energy spectra of 
relic gravitons produced by the pumping action of the gravitational field can very well extend 
for (approximately) twenty five orders of magnitude in frequency . From the physical point 
of view, this observation implies that the energy spectra of relic gravitons can be extremely 
relevant in order to probe the past history of the Universe in a regime which will never be 
directly accessible with observations of electromagnetic backgrounds of cosmological origin 

§• 

In spite of the fact that the nature of the production mechanism is shared by differ- 
ent types of models [|I| , the specific amplitudes of the energy spectra can very well change 
depending upon the behavior of the background evolution. An example in this direction 
are logarithmic energy spectra increasing in frequency Q or even non-monotonic spectra 
originally discussed by Gasperini f5j. 

Different theoretical signals (with different spectral distributions) lead to detector out- 
puts of different amplitudes. Therefore, in order to evaluate the performances of a given 
detector one has to choose the specific functional form of the logarithmic energy spectrum. 
String cosmological models |7j are an interesting theoretical laboratory leading usually to 
sizable theoretical signals in the operating window of wide band interferometers (WBI) |5|. A 
possible detection of these backgrounds would represent an interesting test for cosmological 
models inspired by the low energy string effective action. Other possible choices are repre- 
sented by scale invariant spectra || || or by tilted ("blue" [|ITJ]) spectra whose energetical 



content is typically concentrated at frequencies larger than the mHz ITT] , |T2| . 

The signal induced in the detector output by a stochastic background of gravitational 
radiation is indistinguishable from the intrinsic noise of the detector itself. This implies that, 
unless the amplitude of the signal is very large, the only chance of direct detection of these 
backgrounds lies in the analysis of the correlated fluctuations of the outputs of, at least, two 
detectors affected by independent noisesQ 

The problem of the optimal processing of the detector outputs required for the detection 
of the stochastic background has been considered by various authors [TBI, III] and it was also 



1 We stress that (see, for instance the second paper in Ref. |l^]) the ratio between the minimum signal 
detectable in the cross-correlation and the minimum signal detectable in the case of a single detector is 
~ l/\/T A/ where T represents the duration of the correlation experiment and A/ is the width of the 
frequency band probed by the detectors. This means that for a measurement with A/ = 1 kHz (like in the 
case of the VIRGO interferometer) and T = 1 yr , the minimum signal detectable with a correlation is 10 5 
times smaller than in the case of a single detector. 
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reviewed in Ref. Ul5f . Suppose, indeed, that the signal registered at each detector can be 
written as (we limit ourselves to the case of two detectors (i = 1, 2)) 



Si = hi(t) + m(t) , (1.1) 

where we have indicated with n the intrinsic noise of the detector, and with h the gravita- 
tional strain due to the stochastic background. By assuming that the detector noises are 
stationary and uncorrelated, the ensemble average of their Fourier components satisfies 

«(/wf)> = ^(/-m^G/i), (i-2) 

where S n (\f\) is usually known as the one-sided noise power spectrum and is expressed in 
seconds. Starting to the signals S\ and a correlation "signal" for an observation time T 
can be defined in the following way: 



/•T/2 rT/2 

S= dt dt' si(t) s 2 (f) Q(t - t') (1.3) 

J-T/2 J -T/2 



where Q is a filter function that depends only by t — t' because we assume that n and h 
are both stationary. The optimal choice of Q corresponds to the maximization of the signal- 
to-noise ratio associated to the "signal" S. In this calculation the stochastic background, 
besides to be stationary, is also assumed to be isotropic, unpolarized and Gaussian. Under 
the further assumptions that detector noises are Gaussian, much larger in amplitude than 
the gravitational strain and statistically independent on the strain itself, it can be shown 
14| , |15| that the signal-to-noise ratio in a frequency range (f m , /m) is given by|: 



where Hq is the present value of the Hubble parameter and F is a numerical factor depend- 
ing upon the geometry of the two detectors. In the case of the correlation between two 
interferometers F = 2/5. In Eq. (|1 .4j ) , the performances achievable by the pair of detectors 
are certainly controlled by the noise power spectra (NPS) S^ 1,2 \ However in Eq. (1.4), on 



top of NPS, there are two important quantities. The first one is the theoretical background 
signal defined through the logarithmic energy spectrum (normalized to the critical density 
p c ) and expressed at the present (conformal) time T] [] 

Ogw(/, = - ^7 = nfab) r)o) • (1.5) 
Pc din/ 



2 For a clear discussion about these assumptions see Ref. [jl5| . In this paper are also discussed the 
modifications needed in the case in which most of these assumptions are relaxed. Finally, in order to avoid 



possible confusions we stress that the definition of the SNR is the one discussed in, 1 12 and it is essentially 
the square root of the one discussed in 111 



14 



If 



3 In most of our equations we drop the dependence of spectral quantities upon the present time since all 
the quantities introduced in this paper are evaluated today. 
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As we can see, we have chosen to parametrize the theoretical spectrum through a time- 
independent amplitude (^(770)) and a frequency-dependent part (u(f)). The important 
quantity appearing in eq. ( |1.4|) is the overlap reduction function j(f) 0, |15| which is a 



dimensionless function describing the reduction in the sensitivity of the two detectors (at a 
given frequency /) arising from the fact that the two detectors are not in the same place 
and, in general, not coaligned (for the same location and orientation 7(/) = 1). Since the 
overlap reduction function cuts-off the integrand of Eq. ( |1.4j ) at a frequency which approxi- 
mately corresponds to the inverse separation between the two detectors, it may represent a 
problematic (but controllable) element in the reduction of the sensitivity of a given pair of 
detectors. 

Various ground-based interferometric detectors are presently under construction (GEO 
Tj|, LIGO-LA, LIGO-WA fT7[, TAMA @, VIRGO |TJ). Among them, the pair consisting 



of most homogeneous (from the point of view of the noise performances) detectors with 
minimum separation is given by the two LIGOs (VIRGO and GEO are even closer, but they 
have different performances for what concerns the NPS). However, this separation (~ 3000 
km) is still too large. The overlap reduction function 7(/) for the pair LIGO-LA— LIGO-WA 
encounters its first zero at 64 Hz, falling off (swiftly) at higher frequencies, i.e., right in the 
region where the two LIGOs, at least in their initial version, have better noise performances. 

The motivation of our exercise is very simple. We want to study the sensitivity of a sys- 
tem of two VIRGO-like detectors to stochastic backgrounds of gravitational radiation. Up 



to now, if we exclude the case of Ref . |H| , the correlation between two VIRGO detectors has 
not been seriously explored in the literature. The reason for this lack of studies is that, in 
contrast with the LIGO project where two detectors are simultaneously under construction, 
only one VIRGO detector is currently being built. However, recently, within the European 
gravitational wave community, the possibility of building in Europe an interferometric detec- 



tor of dimensions comparable to VIRGO has received close attention [20]. Therefore, there 
is a chance that in the near future the VIRGO detector, now under construction at Cascina 
(Pisa) in Italy, will be complemented by another interferometer of even better performances 
very close (at a distance d < 1000 km) to it. In this paper we examine in detail the possible 
improvements in the VIRGO sensitivity as a result of direct correlation of two VIRGO-like 
detectors. Technological improvements in the construction of the interferometers can be 
reasonably expected in the next years. Thus the VIRGO detectors will gradually evolve 



towards an advanced configuration [ 20| . For this reason we also examine the possible conse- 
quences of a selective improvements of the noise characteristics of the two detectors on the 
obtained results. In order to make our analysis concrete we will pay particular attention 
to the evaluation of the performances of a pair of VIRGO detectors in the case of string 
cosmological models 0, || . 

In order to evaluate precisely the performances of a pair of VIRGO detectors we will use 
the following logic. First of all we will pick up a given class of theoretical models which look 
particularly promising in view of their spectral properties in the operating window of the 
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WBI. Secondly we will analyze the signal-to-noise ratios for different regions of the parameter 
space of the model. Finally we will implement some selective reduction of the noises and 
we will compare the results with the ones obtained in the cases where the noises are not 
reduced. We will repeat the same procedure for different classes of models. 

The results and the investigations we are reporting can be applied to spectra of arbitrary 
functional form. The only two requirements we assume will be the continuity of the loga- 
rithmic energy spectra (as a function of the present frequency) and of their first derivative. 
We will also give some other examples in this direction. 

The plan of our paper is then the following. In Section II we introduce the basic semi- 
analytical tecnique which allows the evaluation of the SNR for a pair of WBI. In Section III we 
will evaluate the performances of a pair of VIRGO detectors in the case of string cosmological 
models. In Section IV we will show how to implement a selective noise reduction and we will 
investigate the impact of such a reduction in the case of the parameter space of the models 
previously analyzed. Section V contains our final discussion and the basic summary of our 
results. 



2 SNR evaluation 

In the operating window of the VIRGO detectors the theoretical signal will be defined 
through the logarithmic energy spectrum reported in Eq. (|1.5|) . In the present Section 
we shall not make any specific assumption concerning w(f) and our results have general 
applicability. We will only assume that it is a continuous function of the frequency and we 
will also assume that its first derivative is well defined in the operating window of WBI. This 
means that 00(f) can be, in principle, a non-monotonic function. 



2.1 Basic Formalism 



The expected noise power spectrum for the VIRGO detector |T3 is well approximated by 



the analytical fit of Ref . [ 21 1 , namely 
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In order to compute reliably (and beyond naive power counting arguments) the SNR we 
have to specify the overlap reduction function 7(/). The relative location and orientation 
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of the two detectors determines the functional form of 7(/) which has to be gauged in 
such a way that the overlap between the two detectors is maximized (i.e. 7(/) — 1 for 
most of the operating window of the two VIRGO). Moreover, the two interferometers of 
the pair should also be sufficiently far apart in order to decorrelate the local seismic and 
electromagnetic noises. Since the precise location of the second VIRGO detector has not been 
specified so far P0| , we find useful to elaborate about this point by computing the overlap 
reduction functions corresponding to two coaligned VIRGO interferometers with different 
spatial separations. The results of these calculations are reported in Fig. [IQ. Needless to 
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Figure 1: We report the overlap reduction function(s) for the correlation of the VIRGO 
detector presently under construction in Cascina (43.6 N, 10.5 E) with a coaligned interfer- 
ometer whose (corner) station is located at: A) (43.2 N, 10.9 E), d = 58 km (Italy); B) 
(43.6 N, 4.5 E), d = 482.7 km (France); C) (52.3 N, 9.8 E), d = 958.2 km (Germany). The 
third site (C) corresponds to the present location of the GEO detector. Notice that from A 
to C the position of the first zero of 7(/) gets shifted in the infra-red. 

say that these choices are purely theoretical and are only meant to illustrate the effects of 
the distance on the performances of the VIRGO pair. 

The curves labeled with A, B, and C shown in Fig. [I] correspond to different distances 
d between the site of the VIRGO detector (presently under construction in Cascina, near 
Pisa) and the central corner station of a second coaligned VIRGO interferometer. Let us 
now look at the position of the first frequency /, for which j(fi) = for each of the curves. 
We can notice that by increasing d (i.e., going from A to C) the value of /, gets progressively 

4 For illustrative purposes, we assumed that a distance of about 50 km is sufficient to decorrelate local 
seismic and e.m. noises. Such a hypothesis is certainly justified within the spirit of our exercise. 
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shifted towards lower and lower frequencies, linearly with d. This means that, for the specific 
purpose of the detection of a stochastic background of gravitational radiation, the position 
of the first zero of the overlap reduction function cannot be ignored. In general we would 
like fi to be slightly larger than the frequency region where the sensitivity of the pair of wide 
band detectors is maximal. In the explicit examples presented in this paper we will focus 
our attention on the case A. The other two configurations have been the subject of a related 
investigation p2f. 



2.2 SNR and bounds on the graviton spectrum 

By inserting the parametrization ( p..5| ) into Eq. ( |1.4| ) we can write 



SNR 2 = JEo_ -J— J , (2.2) 
5V2V f^S 



where we introduced the (dimension-less) integral 

r^M ^ 7 2 (/' z/)^ 2 (/ z/ 



J2 = I ^ : ^ 



Here the integration variable is v = f/fo, with fo a generic frequency scale within the 
operating window of the interferometer, and the integration domain is restricted to the 
region / m < / < / M (i.e., z/ m < v < u M ). In the following we will choose f = 100 
Hz and, taking into account the frequency behavior of •y(f) (see Fig. [[]), we can assume 
/m = 10 kHz (i.e., z/m = 100). The lower extreme f m is put equal to the frequency fb 
entering Eq. (Q (i.e., v m = 0.02). 

For the chosen values of f and S (see Eq. ( |2.1|) ) one has (H = 100 x h km • s _1 ■ 
Mpc" 1 ): 



. 2 - 4.0 x 10~ 7 ( 1 yr \ 1/2 



K Q ~ - - -jr j SNR . (2.4) 

Since we will often refer to this formula we want to stress its physical meaning. Suppose 
that the functional form of ui(f) is given. Then the numerical value of the integral J can be 
precisely computed and, through Eq. (|2.4|), £1 can be estimated. This quantity, inserted in 
Eq. (|1.5|) , determines for each frequency / the minimum Qqw detectable (for an observation 
time T, with a signal-to-noise ratio SNR) by the correlation of the two detector outputs. 

In the next section, will be compared with two other quantities: Q and f2 max . The 
first is the theoretical value of the normalization of the spectrum, while the second represents 
the largest normalization compatible with the phenomenological bounds applicable to the 
stochastic GW backgrounds. These quantities are of different nature and in order to be more 
precise let us consider an example. 

Suppose, for simplicity, that we are dealing with a logarithmic energy spectrum which 
is a monotonic function of the present frequency. Suppose, moreover, that the spectrum 
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decreases sufficiently fast in the infra-red in order to be compatible both with the pulsar 
timing bound and with the CMB anisotropies bounds. Then the most relevant bound will 
come, effectively, from Big-Bang nucleosynthesis (BBN) p3| [24], |2~5 |. If at BBN too many 



massless particles are present they would cause a faster expansion of the Universe. If the 
expansion would be too fast, then, the correct abundances of the light elements ( 4 He, 3 He, 
Li, D) observed in galaxies could not be reproduced. Thus we should require that the total 
number of massless particles present in the plasma at BBN should not exceed the energy 
density stored, at the same epoch, in radiation. Therefore, in our particular case, we will 
have that fi max is determined by demanding that |23| 



hi J n GW {f,Vo) din/ < 0.2/^(7/0) ^ 5 x 10" 6 , (2.5) 

where Q, 7 (rjo) = 2.6 x 10~ 5 h$ 2 is the present fraction of critical energy density stored in 
radiation. According to our definition, Q max is the maximal normalization of the spectrum 
compatible with the previous inequality, namely, 

ma^ 5 X 10 -6 /"/max 

h 2 n max ~ ' , Z = / co(f) din/. (2.6) 



Z Jf n 

Notice that / ns ~ 1CT 10 Hz is the present value of the frequency corresponding to the horizon 
at the nucleosynthesis time; / max stands for the maximal frequency of the spectrum and it 
depends, in general, upon the specific theoretical model. If the spectrum has different slopes, 
fi max will be determined not only by the nucleosynthesis bound but also by the combined 
action of the CMB anisotropy bound [f| ^ and of the pulsar timing bound |28| . Indeed, we 



know that the very small fractional timing error in the arrival times of the millisecond pulsar's 
pulses implies that Oqw ~ 1CT 8 for a frequency which is roughly comparable with the inverse 
of the observation time along which pulsars have been monitored (i.e., u p ~ f/Z Q b s = 10 -8 
Hz) . Moreover, the observations of the large scale anisotropies in the microwave sky |26|, 



imply that the graviton contribution to the integrated Sachs- Wolfe effect has to be smaller 
than (or at most of the order of) the detected amount of anisotropy. This observation 
implies that f2cw < 6.9 x 10~ n for frequencies ranging between the typical frequency of 
the present horizon and a frequency thirty of forty times larger. In the case of a logarithmic 
energy density with decreasing slope the f2 max will be mainly determined by the Sachs- Wolfe 
bound and it will be the maximal normalization of the spectrum compatible with such a 
bound. 

In general , we will have that f2 < fi max , namely the theoretical normalization of 
the spectrum is bounded, from above, by the maximal normalization compatible with all 
the phenomenological bounds. Therefore, the mismatching between these quantities can 
be interpreted as an effective measure of the theoretical error in the determination of the 
absolute normalization of the spectrum. 

Since 00(f) enters (in a highly non-linear way) into the form of J (as defined in Eq. 
( 2.3 )), the corresponding Q in Eq. fl2.4|) will be different for any (specific) frequency depen- 



dence in u(f). The consequence of this statement is that it is not possible to give a general 



7 



(and simple) relation between the sensitivity at a given frequency, the spectral slope and 
the (generic) theoretical amplitude of the spectrum. However, given the form of the theo- 
retical spectrum, the phenomenological bounds (depending upon the theoretical slope) will 
fix uniquely the theoretical error and the maximal achievable sensitivity. So, if we want to 
evaluate the performances of the VIRGO pair we should pick up a given class of theoretical 
models (characterized by a specific functional form of uj(f)) and compute the corresponding 
sensitivity. The same procedure should then be repeated for other classes of models and, 
only at the end, the respective sensitivities can be compared. 



3 Primordial gravitons versus VIRGO*VIRGO 

We can consider, in principle, logarithmic energy spectra with hypothetical analytical forms 
and arbitrary normalizations. If the logarithmic energy spectrum is either a flat or a de- 
creasing function of the present frequency ||, we can expect, in general, that the theoretical 
signal will be of the order of (but smaller than) 1CT 15 [^7], 29] for present frequencies compa- 
rable with the operating window of the VIRGO pair. This happens because of the combined 
action of the Sachs- Wolfe bound together with the spectral behavior of the infra-red branch 
of the spectrum produced thanks to the matter-radiation transition. Of course this obser- 
vation holds [T2^] for models where the graviton production occurs because of the adiabatic 
variation of the background geometry JTJ, 0] Q 

In order to have large signals falling in the operating window of the VIRGO pair we 
should have deviations from scale invariance for frequencies larger than few mHz. Moreover, 
these deviations should go in the direction of increasing logarithmic energy spectra. This is 



what happens in the case of quintessential inflationary models [11]. In this case, however 



as we discussed in a previous analysis [|T^], the BBN bound put strong constraints on the 



theoretical signal in the operating window of the VIRGO pair. 

Another class of model leading to a large theoretical signal for frequencies between few 
Hz and 10 kHz is represented by string cosmological models ||, [7], Therefore, in order 
to evaluate the performances of the VIRGO pair and in order to implement a procedure of 
selective noise reduction we will use string cosmological spectra. 



3.1 Minimal models of pre-big-bang 

In string cosmology and, more specifically, in the pre-big-bang scenario, the curvature scale 
and the dilaton coupling are both growing in cosmic time. Therefore the graviton spectra 

5 An exception to this assessment is represented by cosmic strings models leading to a flat logarithmic 
energy spectrum for frequencies between fO -12 Hz and 1CP 8 Hz ]|(| plfl . Another possible exception is 
given by the gravitational power radiated by magnetic (and hypermagnetic) ]3^ | knot configurations at the 
electroweak scale Q . 



will be increasing in frequency instead of decreasing as it happens in ordinary inflationary 
models. 

In the context of string cosmological scenarios the Universe starts its evolution in a very 
weakly coupled regime with vanishing curvature and dilaton coupling. After a phase of 
sudden growth of the curvature and of the coupling the corrections to the tree level action 
become important and the Universe enters a true stringy phase which is followed by the 
ordinary radiation dominated phase. It should be stressed that the duration of the stringy 
phase is not precisely known and it could happen that all the physical scales contained within 
our present Hubble radius crossed the horizon during the stringy phase as pointed out in || . 

The maximal amplified frequency of the graviton spectrum is given by 0. § 

/ 10 3XV12 

MVo) ~ 64.8 V^T — GHz (3.1) 



where n r is the effective number of spin degrees of freedom in thermal equilibrium at the 
end of the stringy phase, and g\ = M s /Mp\ where M s and Mp\ are the string and Planck 
masses, respectively. Notice that gi is the value of the dilaton coupling at the end of the 



stringy phase, and is typically of the order of 10~ 2 -j- 1CT 1 [|34|. As we can see from the 
previous equation the dependence upon n r is quite weak. In order to red-shift the maximal 
amplified frequency of the spectrum from the time r\\ (which marks the beginning of the 
radiation dominated evolution) up to the present time we assumed that the cosmological 
evolution prior to t)q and after rji is adiabatic. Minimal models of pre-big-bang are the ones 
where a dilaton dominated phase is followed by a stringy phase which terminates at the 
onset of the radiation dominated evolution. In the context of minimal models, the function 
u>(f) introduced in Eq. ( |1.5| ) can be written as 

3 



"(/) 




2/3-3 



-^4 

2 fs 



/</. = - 



f s < f < h 



(3.2) 



where, 







In z„ 



(3.3) 



In this formula z s = fi/f s and g s are, respectively, the red-shift during the string phase and 
the value of the coupling constant at the end of the dilaton dominated phase. The first of 
the two branches appearing in Eq. ( |3.2| ) is originated by modes leaving the horizon during 
the dilaton dominated phase and re-entering during the radiation dominated phase. The 
second branch is mainly originated by modes leaving the horizon during the stringy phase 
and re-entering always in the radiation dominated phase. The theoretical normalization 



Q th = 2.6 g 2 , (^-) ' 0,(770) 



(3.4) 
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multiplied by uj(f) (as given in Eq. ( |3.2|) ) leads to the theoretical form of the spectrum. 
Notice that n r is of the order of 10 2 4- 10 3 (depending upon the specific string model) and 
it represents a theoretical uncertainty. 

However, as anticipated in the previous section, the theoretical normalization of the 
spectrum should be contrasted with the one saturating the BBN bound (i.e., f2 max ). This 
quantity is obtained by Eq. (|2.6|), where in the case under consideration 



X — Xd + X s 



with 



X d = f fs %u(f), X s = I 1 " %oo{f). (3.5) 

J/ns J Jfs J 



h df 



In the case of minimal pre-big-bang models the analytical expressions of Xd and X s are given 
by 



X/ 
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X, 
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r 2/3-6 „-2/3 



2/3(3-2/3) 2/3-6 2(3 ' 



(3.6) 



In the case of non-minimal models of pre-big-bang the integrals determining the BBN bound 
are instead 



Ji = A(a, z s ) + B(a, z s ) \az s + C(a, z s ) In 2 z s , 



where and z r = f\j f r and 
A(a, z s ) = 



-4 _ -4 



1 + lnz s f 



(3.7) 
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In the intermediate frequency region of the graviton spectra an important bound comes 
from the pulsar timing measurements. Therefore, if one ought to consider rather long stringy 
phases (i.e., large z s ), the BBN constraint should be supplemented by the requirement that 
^gw(10~ 8 Hz) < 1(T 8 |||. We will come back to this point later. 



Following the explicit expression of the function oo(f), Eq. ( |2.4| ) can be re- written as 
follows: 

where, introduced the following notation 

J k = [ Us dv ^ ^2 ln fc z/, k = 0,1,2,3,4 



J±m(3-2f3) ~ / du ^(l), f x v (2),, ^ ; , m-1,2 (3.10) 

C d = l + 2 f- 3 + ilnz/ s , 



one has 



2 s~ 2/3 / 4 3 2 1 1 \ 1/2 

= — 3~ ( Crf^O — 2C d Ji + -C d J2 — -CdJj, + 77^4 



2 a 2 lg 



1/2 



7 - f fi 7 _l 1 ^Z 3 ^ 6 , 4 J_ 4 2/3 ~ 3 I / > I I ■ 

s ~ 7^ + + + 4 ^2£ + 4 ^^3 • l^ 11 ) 

"l \ ^1 ^1 ^1 / 

The previous expressions are general in the sense that they are applicable for a generic 
value of f s . If f m < f s < /m then both J s and Jd give contribution to the sensitivity. If, on 
the other hand f s < f m (i.e., a long stringy phase) the main contribution to the sensitivity 
will come from J s . The integrals appearing in J^ s have to be evaluated numerically. In all 
our calculations we will assume that both VIRGO detectors are characterized by the same 
(rescaled) NPS (reported in Eq. (|2~2|) ). 

The main steps of our calculation are the following. We firstly fix gi and for each pair 
( z s, Qi/ 9s) (within the range of their physical value) we compute Vt (for T = 1 yr and SNR 
= 1), and O max . We then compare these two quantities to the theoretical normalization 
given in Eq. ( p.4| ). If f2 will be larger than f2 (but smaller than fi max ) we will say that 
the theoretical signal will be "visible" by the VIRGO pair. In this way we will identify in 
the plane (z s , gi/g s ) a visibility region according to the sensitivity of the VIRGO pair. The 
theoretical error on the border of this region can be estimated by substituting f2 max to Q th . 

To illustrate this point we consider a specific case. The value of the coupling at the end 
of the stringy phase can be estimated to lie between 0.3 and 0.03 [[34]]. The knowledge of 
gi will not fix uniquely the theoretical spectrum which does also depend on the number of 
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relativistic degrees of freedom at the end of the stringy phase. Therefore, the theoretical 
error in the determination of the absolute normalization of the spectrum could be also viewed 
as the error affecting the determination of n T . In all the plots shown we will take, when not 
otherwise stated, g\ = 1/20 and n r = 10 3 as fiducial values. Different choices of g\ will 
lead to similar results. We will also assume that the overlap reduction function associated 
with the pair is the one reported in the curve A of Fig. [TJ 




2 4 6 B 10 12 14 16 IB 2 4 6 8 tO 12 14 16 II 



log z 5 log z s 

Figure 2: We report the ratios Q m£LX /Q (left) and fl th /fl (right) as a function of Q\j g a and 
logZg (Q is calculated for T = 1 yr and SNR = 1). The lower contour plots show the regions 
where these ratios are greater than 1. The shaded area (bottom right) represents the region 
where the combination of the theoretical parameters is such that the corresponding Q does 
not violate the BBN bound. As we can see the visibility region is reduced. The difference 
between the shaded area in the right plot and the one in the left plot measures the error 
made by assuming as normalization of the spectrum not the theoretical one but the maximal 
one compatible with the BBN. The value z s = 10 8 roughly corresponds to f s ~ Jo- 
in Fig. H (top left) we report the result of our calculation for the ratio between Q max 
and SI as a function of gi/g s and logz s . The contour plot (bottom left) shows the region 
of the plane (logz s , gi/g s ) where this ratio is greater than 1, i.e. the maximal visibility 
region allowed by the BBN bound. In the opposite case, i.e., to /SI < 1, the VIRGO 
pair is sensitive to a region excluded by the BBN. In the right part of Fig. ^ we go one 
step further and we plot the ratio between SI and SI. The shaded area in the contour plot 
(bottom right) is the region of the plane (\ogz s , gi/g s ) where the conditions SI /Si > 1 
and O max /f2 > 1 are simultaneously met. The shaded area in this plot defines the visibility 
region of the VIRGO pair assuming the theoretical normalization of the spectrum. From 
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Fig. ||, by ideally subtracting the shaded area of the left contour plot from the shaded area 
of the right contour plot we obtain an estimate of the theoretical error. The results we 
just presented can be obviously recovered for different values of g\ close to one. However, 
if gi gets too small (and typically below 1/25) the visibility area gets smaller and smaller 
eventually disappearing. 

The visibility regions appearing in Fig. |2| extend from intermediate values of z s (of 
the order of 10 8 ) towards large values of z s (of the order of 10 18 ). Notice that for our 
choice of gi, f s can become as small as 10~ 8 for z s of the order of 10 18 . As we recalled in 
the previous Section, this frequency corresponds to the inverse of the observation time along 
which pulsar signals have been monitored and, therefore, for this frequency, a further "local" 
bound applies to the logarithmic energy spectra of relic gravitons. This bound implies that 
f2Gw(10 _8 Hz) < 10~ 8 . In our examples, the compatibility with the BBN bound implies 
also that the pulsar timing constraint is satisfied. Given our choice for gi we can clearly see 
that the visibility regions depicted in Fig. ^| extend for values of g s which can be as small as 
1/160 (or as small as 1/60 in the case of right part of Fig. ^|). 

3.2 Non-minimal models of pre-big-bang 

In the context of minimal models of pre-big-bang, the end of the stringy phase coincides 
with the onset of the radiation dominated evolution. At the moment of the transition 
to the radiation dominated phase the dilaton sits at its constant value. This means that 
gi ~ 0.03 -T- 0.3 at the beginning of the radiation dominated evolution. As pointed out 
in ||, it is not be impossible to imagine a scenario where the coupling constant is still 
growing while the curvature scale starts decreasing in time. In this type of scenario the 
stringy phase is followed by a phase where the dilaton still increases, or, in other words, the 
coupling constant is rather small at the moment where the curvature starts decreasing so 
that ft < 1. 

After a transient period (whose precise duration will be fixed by the value of gi), we will 
have that the radiation dominated evolution will take place when the value of the coupling 
constant will be of order one (i.e., g r ~ 1). An interesting feature of this speculation is 
that the graviton spectra will not necessarily be monotonic (as the ones considered in the 
previous analysis). We then find interesting to apply our considerations also to this case. 

The function u(f) in the non- minimal model described above is given by || 

• fr < f < fs = ~ 
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• f s < f < h 

"(/) = 

where, in the present case 




2/V3 





2+cr 



i-^jr) ( 3 - 12 ) 



/l ^ 648 ^ I (|) 1/2 '(^"J GHZ ' fr= { 9 fi ] h ' (3,13) 




(Notice that the form of u(f) reported in || differs from our expression only by logarithmic 
correction whose presence is, indeed, not relevant.) The frequency f r corresponds to the onset 
of the radiation dominated evolution. If we adopt a purely phenomenological approach we 
can say that f r has to be bounded (from below) since we want the Universe to be radiation 
dominated not later than the BBN epoch. Thence, we have that f r > / ns . Recalling 
the value of the nucleosynthesis frequency and assuming that g r ~ 1 this condition implies 
gi > 8.2 x 1CT 16 . This simply means that in order not to conflict with the correct abundances 
of the light elements we have to require that the coupling constant should not be too small 
when the curvature starts decreasing. Notice that for frequencies / < f r the spectrum 
evolves as f~ 3 . The ultra-violet branch of the spectrum is mainly originated by modes 
leaving the horizon during the stringy phase and re-entering when the dilaton coupling is 
still increasing. 

Concerning the non-minimal spectra few comments are in order. Owing to the fact 
that gi can be as small as 10~ 15 we have that the highest frequency of the spectrum can 
become substantially smaller than in the minimal case. Moreover, the spectrum might also 
be non-monotonic with a peak at f s . Looking at the analytical form of the spectrum we see 
that this behavior occurs if a > 2. We remind that a parametrizes the spectral slope in 
the phase where the curvature scale decreases but the dilaton coupling is still growing |5|. 
A non-monotonic logarithmic energy spectrum (with a maximum falling in the sensitivity 
region of the VIRGO pair) represents an interesting possibility. 

The results of our calculation for g\ = 10~ 12 , n r = 10 3 , g r = 1, and a > 2 are reported 
in Fig. |3|. As done in the case of minimal spectra we analyse the visibility window in the 
plane of the relevant parameters of the model. As we can see from the left part of Fig. 
|3] the region compatible with the BBN is rather large but it shrinks when we impose the 
theoretical normalization ( right part of Fig. |3|) which is always smaller than the maximal 
normalization allowed by BBN. 

It is interesting to compare directly the three dimensional plots appearing in Fig. [2] 
with the corresponding three dimensional plots of Fig. ^. Notice that the shaded region in 
the case of minimal models corresponds to ratios Q max /Q and Q th /Q which can be 3 or 2, 
respectively. On the other hand the shaded region in the case of Fig. ||| corresponds to ratios 
O /H and f2 /Vt which can be, respectively, as large as 50 or 25. So, in the latter case the 
signal is larger for a smaller region of the parameter space. 
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Figure 3: In order to make clear the comparison with the visibility region of the minimal 
models, we report fi max /f2 (left) and VL /Vt (right) as a function of o and of the logz s in 
the non-minimal scenario. Notice that we took gi = 1CT 12 , n r = 10 3 , and g r = 1. As for 
Fig. 0, the shaded areas in the lower contour plots represent the region where each ratio is 
greater than 1, and, in the case of the right plot, also the BBN is satisfied. 

4 Noise reduction and the visibility region of a VIRGO 
pair 

There are two ways of looking at the calculations reported in this paper. One can look at 
these ideas from a purely theoretical perspective. In this respect we presented a study of the 
sensitivity of a pair of VIRGO detectors to string cosmological gravitons. There is also a 
second way of looking at our exercise. Let us take at face value the results we obtained and 
let us ask in what way we can enlarge the visibility region of the VIRGO pair. In this type 
of approach the specific form of graviton spectrum is not strictly essential. We could use, in 
principle, any motivated theoretical spectrum. As we stressed, we will use string cosmological 
spectra because, on one hand, they are, in our opinion, theoretically motivated and, on the 
other hand, they give us a signal which could be, in principle detected. Of course, there 
are other well motivated spectra (like the ones provided by ordinary inflationary models). 
However, the signal would be, to begin with, quite small. 

In this Section we will then consider the following problem. Given a pair of VIRGO 
detectors, we suppose to be able, by some means, to reduce, in a selective fashion, the 
contribution of a specific noise source to the detectors output. The question we ought to 
address is how the visibility regions will be modified with respect to the case in which the 
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selective noise reduction is not present. We will study the problem for the pair of VIRGO 
detectors considered in the previous Sections, i.e., for identical detectors with NPS given in 
Eq. and characterized by the overlap reduction function of the case A of Fig. |I|. As for 

the theoretical graviton spectrum we will focus our attention on the case of minimal models 
considered in Section III. A, with the same parameters used to produce Fig. ^|. Also here, 
the quantity Q will be computed for T = 1 yr and SNR = 1. 

As shown in Section II the NPS is characterized by three dimension-less numbers S 12i3 , 
and two frequencies f a and Roughly, Si and E 2 control, respectively, the strength of 
the pendulum and pendulum's internal modes noise, whereas S 3 is related to the shot noise 
(see Ref. for an accurate description of the phenomena responsible of these noises). 



Below the frequency /& the NPS is dominated by the seismic noise (assumed to be infinity in 
Eq. (2J.)). The frequency f a is, roughly, where the NPS gets its minimum. The frequency 
behavior of this three contributions and of the total NPS is shown in Fig. |4j. The stochastic 
processes associated with each source of noise are assumed to be Gaussian and stationary. 
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Figure 4: The analytical fit of the rescaled noise power spectrum E n defined in Eq. (|Q| ) 
in the case of the VIRGO detector. With the full (thick) line we denote the total NPS. We 
also report the separated contribution of the three main (Gaussian and stationary) sources 
of noise. 

In the following, without entering in details concerning the actual experimental strategy 
adopted for the noise reduction, we will suppose to be able to reduce each of the coefficients 
Sj by keeping the other fixed. In order to make our notation simpler we define a "reduction 
vector" 

P = (Pi,P2,Pa) , (4.1) 
whose components define the reduction of the pendulum (pi), pendulum's internal modes 
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(p 2 ) and shot (p 3 ) noises with respect to their fiducial values appearing in Eq. (|2.1| ) (corre- 
sponding to the case p = (1, 1, 1)). 

As shown in Fig. ^ the pendulum noise dominates the sensitivity of the detectors in the 
low frequency region, namely below about 40 Hz. In Fig. [| we report the results of our 
calculation for the case p = (0.1, 1, 1). Here the parameters of the theoretical spectrum are 
exactly the same as in Fig. The only change is given by a reduction of the pendulum 
noise. From the comparison between Fig. |^ and Fig. 0, we see that the visibility region in 
the parameter space of our model gets immediately larger especially towards the region of 
small g s . This enlargement is quite interesting especially in terms of Q /Q. 




log z a log z s 

Figure 5: We report the ratios fi max /f2 (left plots), and f2 th /f2 (right plots) in the case in 
which the shot noise and the noise related to the pendulum's internal modes are not reduced, 
whereas the pendulum noise is diminished by a factor of ten with respect to the values quoted 
inEq. (P), i.e.,p = (0.1,1,1). 

In the frequency region between 50 and 500 Hz the performances of the detectors are, 
essentially, limited by the pendulum's internal modes noise. The results obtained for a 
selective reduction of this component are summarized in Fig. where the pendulum and 
shot noises are left unchanged but the internal modes component is reduced by a factor of 
ten. As we can see the visibility region gets larger and the increase in the area is comparable 
with the one obtained by selecting only the pendulum noise. 

Finally, for sake of completeness, we want to discuss the case of the shot noise, i.e., the 
noise characteristic of the detector above 500 Hz. If the shot noise is reduced by one tenth 
(i.e. p3 = 0.1) the visibility region does not increase significantly. This result is consequence 
of the fact that, as shown by Fig. f|, the shot noise contribution to the NPS starts to be 
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Figure 6: We report the result of selective reduction in the case where the noise cause by the 
pendulum's internal modes is reduced by a factor of ten, whereas the pendulum and shot 
contributions are left unchanged, i.e., p = (1,0.1, 1). 

relevant for / ~ 1 kHz, i.e., in a frequency region where the overlap between the detectors 
begins to deteriorate (see Fig. [I]). In Figs. |5| and|^ the thermal noise is reduced by one tenth 
and the increase in the visibility region is, comparatively, larger. This shows that a reduction 
in the shot noise will lead to an effect whose practical relevance is already questionable at 
the level of our analysis. 

Clearly, the simultaneous reduction of both components of the thermal noise leads to a 
substantial increase in the area of the visibility region which gets even larger than the ones 
illustrated in Fig. [| and ||. 



5 Discussion and executive summary 

There are no compelling reasons why one should not consider the appealing theoretical 
possibility of a second VIRGO detector coaligned with the first one. Moreover, recent ex- 



perimental suggestions seem coherently directed towards this goal pDfl . While the location 
of the second detector is still under debate we presented a theoretical analysis of some of the 
scientific opportunities suggested by this proposal. 

We focused our attention on possible cosmological sources of relic gravitons and we 
limited our attention to the case of stochastic and isotropic background produced by the 
adiabatic variation of the backgound geometry. In the framework of these models we can 
certainly argue that in order to have a large signal in the frequency window covered by 
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VIRGO we have to focus our attention on models where the logarithmic energy spectrum 
increases at large frequencies. Alternatively we have to look for models where the logarithmic 
energy spectrum exhibits some bump in the vicinity of the VIRGO operating window. If 
the logarithmic energy spectra are decreasing as a function of the present frequency (as it 
happens in ordinary inflationary models) the large scale (CMB) constraints forbid a large 
signal at high frequencies. In the case of string cosmological models the situation seems more 
rosy and, therefore, we use these models as a theoretical laboratory in order to investigate, 
in a specific model the possible improvements of a possible VIRGO pair. The choice of a 
specific model is, in some sense, mandatory. In fact, owing to the form of the SNR we can 
immediately see that different models lead to different SNR not only because the amplitude 
of the signal differs in different models. Indeed, one can convince himself that two models 
with the same amplitude at 100 Hz but different spectral behaviors between 2 Hz and 10 
kHz lead to different SNR. 

In order to analyze the sensitivity of the VIRGO pair we described a semi-analytical tech- 
nique whose main advantage is to produce the sensitivity of the VIRGO pair to a theoretical 
spectrum of arbitrary slopes and amplitudes. The theoretical error is estimated, in our ap- 
proach, by requiring the compatibility with all the phenomenological bounds applicable to 
the graviton spectra. As an intersting example, we asked what is the sensitivity of a VIRGO 
pair to string cosmological spectra assuming that a second VIRGO detector (coaligned with 
the first one) is built in a european site. By assuming that the second VIRGO detector has 
the same features of the first one we computed the SNR and the related sensitivity achievable 
after one year of observation in the case of string cosmological spectra. 

By using the string cosmological spectra as a theoretical laboratory we then studied 
some possible noise reduction. Our main goal, in this respect, has been to spot what kind of 
stationary and stochastic noise should be reduced in order to increase the visibility region 
of the VIRGO pair in the parameter space of the theoretical models under considerations. 
Our main result is that a selective reduction of each of the three main sources of noise is not 
equivalent. A reduction in the shot noise by a factor of ten does not increase significantly the 
visibility region of the VIRGO pair. A selective reduction of the thermal noise components 
is far more efficient. In particular, we could see that a reduction (of one tenth) of the 
pendulum's internal modes increases the visibility region of four times. The simultaneous 
reduction of the two components of the thermal noise leads to an even more relevant increase. 

The construction of a second VIRGO detector coaligned with the first one and an overall 
reduction of the thermal noise of each detector of the pair leads to what we called "up- 
graded VIRGO" program. The results presented in this paper are obtained in the case 
of a particularly promising class of theoretical models but can be generally applied to any 
logarithmic energy spectrum with similar qualitative results. However, owing to the non- 
linearities present in the evaluation of the SNR it would not be correct assess that they hold, 
quantitatively, without change. We hope that our results and our suggestions may turn out 



to be useful in the actual process of design of the upgraded VIRGO program [20 
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